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Abstract 

Here we report lasing action in limacon-shaped GaAs microdisks with quantum dots (QDs) 
embedded. Although the intracavity ray dynamics is predominantly chaotic, high-Q modes are 
concentrated in the region \ > Xc as a result of wave localization. Strong optical confinement by 
total internal reflection leads to very low lasing threshold. Our measurements show that all the 
lasing modes have output in the same direction, regardless of their wavelengths and intracavity 
mode structures. This universal emission direction is determined by directed phase space flow of 
optical rays in the open chaotic cavity. The divergence angle of output beam is less than 40 degree. 
The unidirectionality proves to be robust against small deviations of the real cavity shape and size 
from the designed values. 
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The microdisk laser is a promising candidate for light sources in integrated photonic 
circuits because of its low lasing threshold and in-plane output . One major drawback 
of a circular microdisk laser is its isotropic output. Cavity shape has been employed as a 
design parameter in optimizing microlaser performance, e.g. deformation from circularity 
to achieve anisotropic output , which is an easy and cost-effective way compared to 

waveguide coupling which requires accurate positioning of a waveguide next to or underneath 
the cavity [f], 7, 8]. However, deformed microdisk lasers produce output beams in multiple 
directions. For many applications, e.g., single-photon sources, unidirectional emission is 
essential. Only a few disk shapes such as spiral and rounded isosceles generate light output 
in a single direction [i], 1(3, 11], but they suffer severe quality (Q) spoiling which leads to very 



high lasing threshold. Employing the characteristic modifications of spatial mode structures 
near avoided resonance crossings was also proposed to achieve both high Q and unidrectional 
emission {l^j. Unfortunately it works only for particular modes and the existence of nearly 
degenerate modes with different far-field patterns smears out the output directionality. 

Utilizing a chaotic microcavity is another possible way [13]. An intuitive expectation for 
modes of a chaotic dielectric cavity is low Q factor and pseudorandom emission pattern, 
because an optical ray undergoing chaotic diffusion in phase space can quickly reach the 
leaky region (where the incident angle x a ^ cavity boundary is less than the critical angle 
Xc for total internal reflection) and escape from the cavity by refraction. However, wave lo- 
calization along; unstable periodic orbits (UPOs) creates scarred resonances, which can have 
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high quality factors if the UPOs are located beyond the leaky region in phase space [14 1151]. 
Emission directionalities of these high-Q scarred resonances are determined by the unstable 
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bids of UPOs, which guide optical rays out of the cavity in certain preferred directions 



171 ]. Hence, an open chaotic microcavity can produce directional emission via phase 
space flow along the unstable manifolds [3]. Moreover, many high-Q scarred resonances 
may have similar far-field patterns if their corresponding UPOs have closely-nested unstable 
manifolds j^l, []j]] . The chaotic microcavity lasers that have been fabricated so far have out- 
put beams in multiple directions [2(1 \2l\ . Wiersig and Hentschel recently showed in their 
calculations that a limacon cavity can support high-Q modes with unidirectional emission. 
Here we report lasing in limacon-shaped GaAs microdisks with InAs quantum dots (QDs) 
embedded. Although the intracavity ray dynamics is predominantly chaotic, the high-Q 
modes are concentrated in the region x > Xc as a result of wave localization. Strong optical 
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confinement by total internal reflection leads to very low lasing threshold, which is compa- 
rable to that of a circular disk of similar size. Our measurements show that all the lasing 
modes have output in the same direction, regardless of their wavelengths and intracavity 
mode structures. This universal emission direction is determined by directed phase space 
flow of optical rays in the open chaotic cavity. The divergence angle of output beam is less 
than 40 degree. The unidirectionality proves to be robust against small deviations of the 
real cavity shape and size from the designed values. 

The sample is grown on GaAs substrate by molecular beam epitaxy. The layer structure 
consists of 1000 nm Alo.6sGao.32As and 265 nm GaAs. Inside the GaAs layer there are 
six monolayers of InAs QDs equally spaced by 25 nm GaAs barriers. Limacon-shaped 
microdisks are fabricated with photolithography and two steps of wet etching. The first step 
is a nonselective etching with HBr, followed by a HF-based selective etching to undercut 
Alo.68Gao.32As. Figure 1 shows the top-view and side-view scanning electron microscope 
(SEM) images of a microdisk. The disk shape is slightly deviated from the limacon of 
Pascal, and it can be fitted as p = R(l + ecos#) — rf, where R = 2.28 /xm, e = 0.43, and d 
= 0.1 /im. The side-view SEM image illustrates smooth sidewall of GaAs disk. The 1000 
nm long Alo.6sGao.32As pedestal separates the GaAs disk from the substrate. To minimize 
its effect on lasing modes in the disk, the pedestal is etched to have a top lateral dimension 
of 620 nm . 

In the lasing experiment, the microdisk is cooled to 10K, and optically pumped by a 
mode-locked Ti:Sapphire laser at wavelength 800 nm. A long-working-distance objective 
lens focuses the pump beam to a single disk from the top. A small amount of emission 
from the microdisk is scattered by the disk boundary to the top. It is collected by the same 
objective lens and directed to a spectrometer. Figure 2(a) shows a spectrum of emission from 
the disk shown in Fig. 1 at the incident pump power of 18.1 /xW. It features several narrow 
peaks. Figure 2(b) plots the intensity and linewidth of one peak at A = 923 nm as a function 
of incident pump power P im When Pi exceeds 2.5 /iW, the peak intensity exhibits a rapid 
increase with pumping, meanwhile the linewidth drops quickly to 0.16 nm. This threshold 
behavior indicates the onset of lasing oscillation. The inhomogeneously broadened gain 
spectrum of InAs quantum dots allows lasing in multiple modes which are well separated in 
wavelength. The pump power required to reach the lasing threshold is comparable to that 
of a circular microdisk with the same area [22]. The low lasing threshold suggests efficient 
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confinement of light in the limacon cavity. 

To measure the far-field pattern of laser emission from a limacon cavity, we fabricate a 
ring structure around each microdisk. The ring is centered at the disk and has a radius of 
34 /im, which is much larger than the disk radius R. The in-plane laser emission from the 
disk propagates to the ring and is scattered out of the plane. The scattered light pattern 
is imaged by the objective lens to a CCD camera. Figure 3(a) shows an image taken at 
Pi = 18.1 /iW. To obtain a clear image of weak scattered light along the ring, we use long 
integration time for the CCD camera, which causes an over-exposure of the disk itself. Since 
the ring radius exceeds 4i? 2 /A, the scattered light intensity along the ring reflects the far-field 
emission pattern of the microdisk. It shows the laser output from the limacon disk is mostly 
in one direction. The red line in figure 3(b) plots the scattered light intensity as a function 
of polar angle 9. The unidirectional emission is centered around 9 = with a width of 40°. 
The emission spectrum taken simultaneously with the image [Fig. 2(a)] reveals multi-mode 
lasing. Hence, the unidirectional emission shown in Fig. 3(a) comes from all lasing modes. 
To find out the far-field pattern of individual lasing modes, we place a narrow bandpass 
filter in front of the CCD camera to select only one mode. By tuning the filter's wavelength, 
we image different lasing modes. The far-field patterns of two lasing modes at 909 nm and 
923 nm are also plotted in Fig. 3(b). They are similar to the far-field pattern of total laser 
emission except for a small variation in the angular distribution of output intensity. Our 
data confirm that all the lasing modes have output beams in the same direction with similar 
divergence angle. 

To identify the lasing modes, we have performed numerical simulations of actual mi- 
crodisks that are measured. Since its radius is much larger than its thickness, the microdisk 
can be regarded as a two-dimensional cavity with an effective index of refraction n e ff. For 
transverse electric (TE) polarization (the electric field parallel to the disk surface), we cal- 
culate n e ff = 3.13. Experimentally the lasing modes are TE-polarized, and correspond to 
high-Q modes of frequencies within the gain spectrum. After obtaining the exact shape of 
the microdisk from the SEM image, we compute the high-Q TE resonances by solving the 
Maxwell's equations with the finite-difference time-domain (FDTD) method. Due to limited 
accuracy in determining the refractive index of GaAs at low temperature, the frequencies 
of high-Q modes do not match exactly those of lasing modes. Figures 4(a) and (b) plot the 
spatial profiles of two high-Q modes in the microdisk shown in Fig. 1. The mode at A = 928 
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nm has Q = 94000, and the one at 945 nm has Q = 15000. Figures 4 (c) and (d) are the 
Husimi phase-space projections of these two modes. It is evident from Fig. 4(d) that the 
mode at 945 nm is localized on an unstable periodic orbit with four bounces from the cavity 
boundary. The locations of the bounces are marked by green dots in Fig. 4(d). The orbit 
itself is drawn in the inset of Fig. 4(d). Hence, the mode at 945 nm is a scar mode. From its 
spatial profile and Husimi projection, the mode at 928 nm is more like a whispering-gallery 
(WG) mode than a scar mode. However, the whispering-gallery trajectories in a limacon 
cavity with e = 0.43 are confined to a tiny region | sinxl > 0.99 [131, while the mode at 928 
nm is concentrated mostly between | sinx| = 0.74 and | sinxl = 0.93. Strictly speaking, the 
mode at 928 nm is not a WG mode. No stable islands exist below the WG region and above 
region. Hence, we believe the mode at 928 nm is formed by dynamical localization 
Namely, the chaotic diffusion in phase space is suppressed by destructive inter- 
ference, which leads to the formation of WG-like mode. As shown in Figs. 4 (c) and (d), 
both the scar mode at 945 nm and the WG-like mode at 928 nm are localized well above 
the leaky region (sinx < l/n e ff). Their exponentially small intensities in the leaky region 
explain their high Q factors. Another consequence of wave localization in the region \ > Xc 
is that the modes have very low field intensities in the disk center. Thus degradation of 
their Q factors by the Alo.6sGao.32As pedestal is negligible. We have calculated many modes 
and confirmed that all the high-Q modes are concentrated well above the leaky region due 
to wave localization. 

Figures 4(e) and (f) plot the angular distributions of far-field intensities for the two 
high-Q modes of A = 928 nm and 945 nm. Although their intracavity mode structures 
are different, the far-field patterns are similar. Their outputs are predominantly in a single 
direction = 0. Our numerical simulations verify that all high-Q modes have outputs in the 
same direction. To understand the universal directionality of high-Q modes, we take a close 
inspection of their Husimi projections in the leaky region. In Figs. 4(c) and (d), the color- 
enhanced Husimi projections below the critical line sinx = l/n e ff reveal that the escape 
routes from the cavity are nearly identical for the scar mode and the WG-like mode. Wiersig 
and Hentschel showed in |l3|] that light escape for high-Q scar modes in the limacon cavity is 
along the unstable manifolds. The UPOs above the leaky region have closely-nested unstable 
manifolds, which cross the critical line at almost the same point. Hence, the output beams 
have identical direction. Our calculations show that the high-Q modes are not necessarily 
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scar modes, they can be WG-like modes formed by dynamical localization. Although they do 
not correspond to specific UPOs, the WG-like modes have similar unidirectional output as 
the scar modes. To confirm the output directionality is determined by chaotic ray dynamics, 
we have performed classical ray tracing in our microdisk whose shape is slightly deviated 
from the limacon. The initial rays, with identical amplitudes, are uniformly distributed in 
the phase space above the critical line. At each reflection from the boundary, the amplitude 
of a ray is reduced according to the Fresnel law. Tracing of one ray is stopped after its 
amplitude falls below a certain value. The accumulated amplitude distribution in the leaky 
region, shown in Fig. 4(g), resembles the intensity distribution of the WG-like mode in 
Fig. 4(c) and the scar mode in Figs. 4 (d). The good agreement between the ray tracing 
result and the Husimi projection of high-Q modes confirms that the universal directionality 
of output beam results from the classical ray dynamics. 

In conclusion, we have achieved simultaneously low lasing threshold and unidirectional 
emission in a chaotic microcavity. The cavity shape is defined by the limacon of Pascal, and 
the cavity size is smaller than 5 /xm. We observe multi-mode lasing thanks to the broad 
gain spectrum of InAs QDs. Wave localization above the leaky region in phase space en- 
sures strong optical confinement and low lasing threshold. All the lasing modes have output 
beams in the same direction with a divergence angle less than 40°. The emission direction- 
ality is determined by chaotic ray dynamics in the open cavity, and is robust against small 
deviation of cavity shape from limacon. Such tolerance allows fabrication with standard 
photolithography and wet chemical etching, that would facilitate parallel mass production. 
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Figure Captions 

Fig. 1: Scanning electron microscope images of a GaAs microdisk on top of a 
Alo.68Gao.32As pedestal, (a) Top view, (b) side view. The GaAs disk is 265 nm thick 
and contains six layers of In As QDs. 

Fig. 2: (a) Emission spectrum taken at the incident pump power of 18.1 /iW. It features 
several narrow peaks, (b) Intensity (black dot) and spectral width (red cross) of one peak 
at A = 923 nm as a function of incident pump power P im A threshold is clearly seen at Pi = 
2.5 /xW where the slope of peak intensity (marked by black line) changes and the spectral 
width drops quickly. 

Fig. 3: Measured directionality of laser emission from the microdisk shown in Fig. 1. 
(a) An image of laser emission from the microdisk scattered by the ring surrounding the 
disk. The ring radius is much larger than the disk radius. The incident pump power is 18.1 
/xW. The emission spectrum taken simultaneously with the image is shown in Fig. 2(a). (b) 
Far-field angular distributions of emission intensities of all lasing modes (red solid line), as 
well as of two lasing modes at A = 909 nm (blue dashed line) and 923 nm (black shot dash 
line), respectively. All lasing modes have output beams in the same direction (9 is the angle 
in limacon fitting formula) with a divergence angle of 40°. 

Fig. 4: Numerical simulation results of high-Q modes in the microdisk of shape equal 
to that in Fig. 1. Calculated spatial intensity distributions (a, b), Husimi phase space 
projections (c, d), and far-field emission patterns (e, f) of two high-Q modes at A = 928 
nm and 945 nm, respectively. The horizontal axis of (c) and (d) represents the length along 
the cavity boundary from the point 9 = normalized by the cavity perimeter and the 
vertical axis corresponds to sinx, where x is the incident angle at the cavity boundary. 
The spatial pattern and Husimi projection of the mode at 928 nm are similar to those of 
a whispering-gallery mode. The mode at 945 nm is localized on an unstable periodic orbit 
with four bounces from the cavity boundary. The bounces are marked by the green dots in 
(d) and the orbit itself is drawn in the inset of (d). In (c) and (d) the intensities of Husimi 
projections in the leaky region sinx < V n e// are enhanced to illustrate the escape route of 
light from the cavity, (g) is the distribution of optical rays in the leaky region obtained by 
classical ray tracing. It agrees well with the Husimi projections of high-Q modes, illustrating 
the universal output directionality results from chaotic ray dynamics in the open cavity. 
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Fig. 3 Song et al. 
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Fig. 4 Song et al. 
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